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ABSTRACT 



> 

o 
en 

r~^ ] We study the Minimal Messenger Model, a minimal version of Gauge Mediated Su- 

0^ I persymmetry Breaking models. Boundary conditions equal to zero for trilinear and 

r~| ' bilinear soft parameters at the messenger scale make this model free from the super- 

Mh, symmetric CP problem and extremely predictive. These boundary conditions and the 

Q^l vicinity of the messenger scale to the electroweak one, requires a careful implementa- 

tion of the mechanism of radiative breaking of SU{2) x U{\). We assess the importance 
of considering the complete one-loop effective potential and of including a set of log- 
arithmic two-loop corrections to the B parameter for the correct determination of 
/\ ' the electroweak minimum. We analyze the resulting low-energy spectrum and give 

j^ . predictions of interest for future experimental searches. 



1. Introduction 

Models with low-energy breaking of supersymmetry, communicated by gauge interactions to 
the observable sector, have recently drawn considerable attention as interesting alternatives to 
models in which this communication is mediated by gravity [1[. The minimal realization of the 
latter at the electroweak scale is the well known Minimal Supersymmetric Standard Model (MSSM). 

The former, the so-called gauge mediated supersymmetry breaking (GMSB) models have, in- 
deed, several attractive features. Most important of these is the fact that, as gauge interactions 
are flavor blind, squark and slepton masses are universal. Given the low value of the supersymme- 
try breaking scale, this universality is hardly broken by the evolution of mass parameters to the 
electroweak scale through Renormalization Group Equations (RGE). (Universality means, in this 
context, that scalar masses are only functions of gauge quantum numbers, and that A-terms are 
small or proportional to fermion yukawa couplings.) Moreover, they can be more predictive than 
the MSSM, having a smaller number of free parameters, and may, at least in a minimal variant Q, 
provide a solution to the supersymmetric CP problem ^. 

In this minimal version, dubbed the Minimal Messenger Model (MMM), trilinear and bilinear 
soft parameters vanish at the messenger scale X. Moreover, if X is of 0{A), with A the ratio of 
the messenger F-term over X, after the radiative breaking of the electroweak gauge symmetry is 
implemented, this model turns out to be practically a one-parameter model, A. 

Although the MMM was already considered in ||2|,Q,|5| , the mechanism of breaking of the elec- 
troweak gauge group was not always correctly implemented. Furthermore, not all experimental 
constraints on the model were always included. The aim of this paper is to present a comprehen- 
sive analysis addressing these issues. After a definition of the model in Sect. 2, we will impose the 
breaking of SU{2) x U{\) through minimization of the RGE improved tree-level Higgs potential 
(Sect. 3.1). In Sect. 3.2, we will demonstrate the importance of considering the full one-loop cor- 
rected effective potential for the determination of the electroweak minimum. We will also argue on 
the need to include a set of additional two-loop corrections to the parameter B of same size than 
those induced by the one-loop effective potential. Finally, we will discuss the viability of the model 
for different values of A. In the last section, Sect. 4, we will verify which regions of A survive the 
imposing of experimental bounds coming from direct searches of supersymmetric particles and the 
indirect constraint due to the measurement of 5 ^ S7 and list the main predictions of this model. 



2. The model 

The messenger sector of the MMM consists of only one pair of chiral superfields, $, ^, which 
transforms as a vectorlike representation of the electroweak gauge group. In particular, having in 
mind an embedding of this model in an underlying Grand Unified Theory (GUT), $, $ are chosen 
to be in a 5, 5 representation of SU{5) D SU{3)c x SU{2)l x U{1)y- They couple at the tree-level 
to a singlet S (W = \S^^). For simplicity, we neglect the difference in the evolution of A for the 
components of the two 5-plets with different SU{3)c and SU{2)l quantum numbers. 

The scalar component of S acquires a vacuum expectation value (VEV ) giving therefore a 
supersymmetric mass A <S> to fermion and scalar components of both messengers. The auxiliary 
component of S acquires a VEV as well, and gives rise to a supersymmetry-violating mass term, 
A-F5, which mixes the scalar components of ^ and $. The information of supersymmetry breaking 



is then transmitted to the visible sector via gauge interactions. A soft mass is generated for the 
i-th gaugino (i = 1,2,3) at the one-loop level, with fermion and scalar components of <1> and •§ 
exchanged as virtual particles. This is: 

Mr (A <S>) = ^a,{\ <S>) ^ g{x) ; x ^ -^-^ , (1) 

where a GUT normalization for the U{1)y coupling was chosen: ai = (5/3)ay = (5/3)a/cos^ ^vi^. 
If we indicate with A the ratio Fs/ <S> and with X the supersymmetric mass A <S>, we can 
rewrite (|l]) in a more compact form as: 

M,{X) = a,{X)Ag{x); a,{X) ^ ^^ . (2) 

A tilde on gauge and yukawa couplings indicates hereafter a division over Air. Natural values for x 
are of 0(1) (i.e. (A <S>)^ ^ A-F5), but at 2; = 1 massless scalar messengers appear in the theory. 
We choose in the following x = 1/2 and we shall comment on modifications obtained in the limit 
x — > 1. In the range 0—1 the function g{x), given in ||6|, is monotonically increasing and has values 
g{0) = 1, g{l) = 1.386. 

Communication of supersymmetry breaking is passed to the scalars of the observable sector at 
the two-loop level with scalar and fermion components of the messenger fields, gauge bosons and 
gauginos exchanged as virtual particles. The masses obtained are: 

m^iX) = 2A2 I^Cs &liX) + C2 aliX) + ^Y^ a?(X)| f{x) , (3) 

where C3 = 4/3,0 for triplets and singlets of SU{3)c, C2 = 3/4,0 for doublets and singlets of 
SU{2)l; Y = Q — T3 is the hypercharge and i runs over all scalars present in the theory. The 
function f{x), derived in Q and [^, is almost always a flat function equal to 1, except for a sharp 
drop at X = 1, where it has the value 0.7. 

Finally, trilinear and bilinear couplings in the soft scalar potential f\ 

Vsoft ^ -{AuKy^HuQiU^ + {AdhdY^ HdQ^D^ + {AeKY^HdU^ - Bfi HdH^ (4) 

(where Q, U^, D^ etc., as well as -fT^ and H^ indicate here the scalar components of the corre- 
sponding superfields) vanish at the messenger scale: 

A:^{X) = A'iiX) = j^i{X) = ; B{X) = . (5) 

The last relation in (^) is the identifying property of this model. (For a discussion on how such a 
boundary condition can be theoretically implemented, see P|.) Phenomenologically, it renders the 
model very predictive; technically, it makes the search of the correct electroweak minimum rather 
complicated. 

Low-scale (M^) inputs of our analysis are: as = 0.120, 02 = 0.0335, ai = 0.0168, corre- 
sponding to a~^ = 127.9 and sin^ 6w = 0.2316, and loosely compatible with a gauge couplings 
unification [Q. As running fermion masses we use mt{Mz) = 171 GeV, mb{Mz) = 3.0 GeV and 
rririMz) = 1.75 GeV. 

*This relation, the explicit form of the superpotential W = hl^HdLiE] + h''^HdQiDj - K^HuQiU] - fi'^HaHu, 
and the multiplication rule HdQ = ^ijHdiQj\ £12 = — £21 = 1, define the sign-conventions used in this analysis 



3. Radiative Breaking of SU{2)l x U{1)y 

3.1. RGE improved tree-level Higgs potential 

We evolve all the initial parameters (y), (y), and dsj), which depend only on A, from X (= 2A) 
to a decoupling scale Qq. Henceforth, we shall refer to X as to the "high-scale", as compared to 
the "low-scale" ~ Mz- The evolution is performed using the one-loop MSSM RGE. For reference, 
we report in Appendix A these equations and the two-loop equation for B, in the approximation of 
Kobayashi-Maskawa matrix if ~ 1. Those actually used for this analysis, with all intergenerational 
mixing terms needed for the calculation of 6 — > S7, can be found in y]. 

Given the much more modest evolution of the masses of weakly interacting sparticles, and the 
fact that they are much lighter than the strongly interacting ones, we take Qq to be the geometrical 
mean of the high-scale values of m~ and m~. These depend on the decoupling scale itself (see 



(|3|)) since the high-scale gauge couplings are to be obtained from our inputs at Mz through a 
running of the Standard Model (SM) RGE from Mz to Qq and the MSSM RGE from Qq to X. A 
simple iteration allows to find this scale rather quickly. For A ~ 100 TeV, Qq is typically ~ 1 TeV. 
All squark masses cluster around this value. Corrections for the inadequacy of this scale for the 
weak mass parameters involved in the breaking mechanism will be provided by the inclusion of 
the one-loop corrections to the scalar potential, leaving therefore an overall scale ambiguity of the 
next order. Unless a different value is explicitly mentioned, the choice of Qq specified above is that 
made throughout this paper. It will appear obvious later on why this is indeed a good choice. 

For the high-scale yukawa couplings required as inputs of the MSSM RGE, we need the value of 
tan/? . This parameter, together with /i, is obtained by imposing that the electroweak minimum is a 
minimum of the neutral Higgs potential, and, indeed, the deepest one. We postpone the discussion 
of the latter point to a later moment. As for the former, we require that the first derivatives with 
respect to the neutral higgses H^, H^ of the RGE improved tree-level potential 



Vo (Qo) = mL |^.f + /^k |^°f - {B^^) (iiX + h.c) + ^^^ {\hT - l^„T)' (6) 



8 
( /i|^. = mjj, + fJ-'^ {i = n, d)) are equal to zero when it is {H^) = Vu and {H^) = va'- 

(nil^+^,^ + l{g^+g^){vj-vl))2vd-{Bf,)2v^ = (7) 

[ml^+^^-l{g^+gl.){vj-vl))2vu-{Bfi)2vd = 0, (8) 



(' 



and that the obtained solution is a minimum. This minimization condition has to be imposed at 
the low-energy scale Mz- Gauge couplings are indeed evolved down to Mz, whereas, as said before, 
the evolution of mass parameters is stopped at Qq: hence, the Qo~dependence of the potential Vq- 
The equations (^ are more often cast in the form: 



'mjr - mi, tan^ /? 1 o 2Bu 

T^ - -^Iz ; sin 2/? = -^ -^ ^ 

tan2 [3-1 2 ^ ' '^ m]j + m]j + 2^i^ 



l^^= "^''\ 7^T" ->4; sin2/?= , /T ,. . . (9) 



Self-consistent solutions of this problem are obtained through a numerical iteration. 

It is a priori not obvious that the previous equations can yield physically meaningful results, 
due to the little freedom which the model allows and the smallness of the logarithms resummed 
by the RGE, i.e. of the period of integration ixQo = 21n(X/(5o)- The peculiar interplay between 
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Figure 1: Gaugino masses as a function of A, solid lines. The dashed lines indicates the corresponding 
high-energy masses Mi{X). On the right frame are shown the ratios rij and Zi defined in the text. 



strongly and weekly interacting sectors in this model, as we shall see, makes possible the breaking 
of the electroweak gauge group; maintains B small at low-energy, therefore inducing large values 
of tan/3 , and has consequences for the viability of the model itself. 

The evolution of supersymmetric parameters for which yukawa couplings do not play any role 
is independent of the breaking mechanism. Low-scale gaugino masses are: 



M, = M,{X) Zi = ASiiQo) g{x) 



Zi 



aijQo) 

aiiX) 



{l-ai{Qo)bitxQo] 



(10) 



where it is: ai(Qo) = 5j/(l — Oj bf^^ tq^z) and tq^z = 21n (Qq/Mz)- (If no scale is specified, it is 
understood that the relevant variables are low-energy variables (~ Mz), with evolution frozen at 
Qo for the massive ones.) For A = 100 TeV, the period of integration txQo i^ 10-6 to be compared 
with the value ~ 60 in the MSSM for the same decoupling mass and the high-scale X coinciding 
with Mgut ~ 3 X 10^^. The values of the three gauge couplings at X and Qo as obtained for 
A = 100 TeV, starting from our input values at Mz, are: 



ai{X) 
ai(Qo) 



0.0191; 
0.0173; 



a2{X) 
a2{Qo) 



The coefficients Zi are then: 0.90,0.97,1.23 



= 0.0331 ; 
= 0.0322; 

1,2,3). 



a3(X) =0.0739, 
Q3(Qo) = 0.0910. 



ill) 



Gaugino masses, evolution coefficients Zi and ratios rij = Mi/Mj = ai{Qo)/aj{Qo) are shown 
in Fig. (|l]) as function of A. Although we analyze this model for all A's from 15 to 150 TeV, 
we show in this figure only the region from A ~ 62 TeV, for reasons which will become clear 
after the inclusion of one-loop corrections to Vq. For lower A's, the values of Zi hardly deviate 
from those given above. Almost unchanged is also ri2, whereas r32 reaches 3.3 at A ~ 20 TeV. 
Again, the values of Zi are to be compared to those obtained in the MSSM for the same scale Qq: 
Zf^ = 0.44,0.84,2.35(z = 1,2,3). For each value of A, the gaugino sector in the MMM can be 
identified with that in the MSSM for M = Kaig{x)/Zf^ . 



The first two generations of scalar masses, whose boundary conditions can be re-expressed in 
terms of low-scale gaugino masses (|3|) as: 

evolve according to ( |A5| ) in Appendix A and get corrected by quantities ^m^ = rrii^ — mf{X): 

For our choice of x, it is f{x)/{g{x)'^) ~ 1. The corrections (|l^) amount to ~ 15% in the squark 
sector and at most 5% in the case of sleptons, modestly enlarging the already wide gap existing 
between squark and slepton spectra. Dropping the indices of first and second generation, we have, 
for A = 100 TeV: 

m~~16.8M|; m| ~ m| ~ 15.4M|; m~ ~ 1.64M|; m| ~ 0.40M| . (14) 

This large gap sharply distinguish the MMM from the MSSM, which in general predicts for the 
first two generations of squarks and sleptons the familiar relations : 

m~ ~ m^+6.4M|; m~ ~ m~ ~ m^ + 5.9M|; m~ ~ m^+0.68M|; m~ ~ m^+0.22M|. 

Q U D L E 

Therefore, whereas the gaugino sector in the two models can be identified through a specific choice 
of M for each value of A, an identification of the scalar mass parameters would require very different 
values of m? for squarks and sleptons (IO.4M2 ^"^^ O.OGMI for the SU{2)l doublets; Q.hM2 and 
0.18M| for the singlets). 

The interplay between such a light "weak" sector and the heavy "strong" one is such to turn 
m?fj to large negative values as in the MSSM, making therefore possible the radiative breaking of 
SU{2)l X U{1)y- If, for the purpose of illustration, we keep only the first of the logarithms to be 
resummed when solving ([A7|), we get: 



m 



l^ = mliX) - 3at{X) (m|(X) + m|(X) + m|(X)) txQ, (15) 



with at{X) = hf{X)/{ATT) and ht{X) of 0(1). The factor ~ 6 lost in txQo with respect to the 

MSSM is compensated in (^) by the heaviness of the squark spectrum. A similar effect si observed 
2 



in the evolution of m? 



Cancellations between weak and a higher-loop strong terms appear also in the determination 
of the low-energy value of B. Starting from the boundary condition (|5[), a value different from zero 
is generated at the one-loop level through gaugino mediated loops (see first diagram in Fig. ^) . The 
leading contributions are resummed in a series J2n ^^ni^^i-^) txQo)^ where a{X) indicates generically 
a2iX) or ai{X). In a closed form, it is: 

--{-(f)(^)4G)(t)(^)}^ 

Values of Ai different from zero are similarly obtained from one-loop diagrams as the second 
one in Fig. ||. Leading logarithms arising from such diagrams, after resummation, give: 

MsX /l-ZgX fM2\ fl-Z2\ , 3^y /Mi\ /I-Z1 



■^--^^^i^Jl^j+^H^Jl^Jn^-lzfJl^)^ <") 
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Figure 2: One-loop diagrams contributing to B^ and At- Here and in the following figure, the cross 
indicates the parameter //, the dots one or two powers of yukawa couplings. 
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Figure 3: Two-loop diagrams generating B^ and At- Both sets of diagrams contribute to the corresponding 
one-loop and two-loop RGE. 



The coefficients CY are given in Appendix A and C3, C2 are equal to 4/3,3/4 for At and A^, and 
0,3/4 for At-. When only gaugino-mediated loops are considered, it is then, At ^ A^, ^ —M2 and 
B r^ Ar O.IM2, roughly for any A ^^ 70TeV. 

Thus, the leading contributions to B proportional to yukawa couplings, i.e. "Aj-induced" , 
are first obtained at the two-loop level (see upper diagram in Fig. ^) and give rise to the series 
Yl,n Cn+i5t(-'^)(S3(X))" {txQoY^^ ■ At the two-Ioop level are also obtained the leading logarithmic 
contributions to Ai proportional to yukawa couplings (see lower diagrams in Fig. H for At). Nev- 
ertheless, while these diagrams produce indeed numerically small corrections to the values of Ai 
previously obtained, the same is not true in the case of B, for which one has: 






Ms 
M2 



Z2 
Zs 



a2iX) 



{at{X)txQo); 



k: 1. 



(18) 



Large cancellations between weak and higher-loop-order strong terms contributing to B take place. 
They induce a flip of sign for this parameter, which, then, turns out to be small and positive. The 
positivity of B forces also fi to be positive; its smallness, relative to the heavy scalar spectrum, 
pushes tan/3 to large values. 



At the electroweak minimum, after minimization of Vq, one obtains tan/3 = 46.4 and the 
following values of yukawa couplings at the three relevant scales X, Qq and Mz'- 



ht{X) = 0.8130 
/ifc(X) = 0.6324 
hr{X) = 0.4812 



ht{Qo) = 0.8755 ; ht{Mz) = 0.9799 , 

/i6(Qo) = 0.7045 ; hb{Mz) = 0.7984 , (19) 

/i^(Qo) = 0.4716 ; hr{Mz) = 0.4657 . 



The Higgs mass parameters m?^ and m^j are: m^ ~ — (1.50M2)^ — — (0.53M3)^, m?^ ~ 
— (0.97M2)'^ and /i, given practically by —m?H i^^^ difference between the two numbers deter- 
mining Mz), \s 11 c^ I.48M2. The values of At and Ai, do not deviate much from those obtained 
according to (|l^: At ^ A^ = -O.9M2. Ar is now -0.055M2, and B only ~ O.O2M2. 

The results for tan/3 -/u obtained from the minimization of Vq, as a function of A are shown 
in Fig. 1^. Some portions of these lines are dotted to indicate that the corresponding points of 
parameter space are non-physical. Negative squared masses are obtained in the spectra relative to 
these points, for scalar other than H^ and Hd. 

The lightest low-energy soft parameter is certainly m~. For small values of A, however, the 

E 

inclusion of D-terms renders sneutrinos lighter than charged sleptons. In the leptonic sector, D- 
terms are: 

D^=^cos2/3M|; D~^ = -D~{1 -2 sir? dw); D~^ = -D~2 sm^ Ow , (20) 

and numerically they amount to D~ ~ — (64GeV)'^, D~ ~ D^ ~ (45GeV)^ for all values of 
tan/? obtained here through minimization of dm. Moreover, if A is small enough, m^ may become 
negative. By using (10)-(|lJ) (i.e. neglecting, for a rough estimate, third generation mass effects on 



the evolution of soft parameters) and our gauge couplings inputs, it is easy to see that the condition 
ml > \D~\, i.e. m^ > 0, is fulfilled for A > 7.3 | cos2/3|^/2^z/a2 ^ 20TeV. This value is not too 
dissimilar from that obtained with an exact calculation, i.e. the value of A in Fig. ^ where the 
initial dotted intervals of the two curves for fj, and tan/3 turn into solid ones. 

In this same lowest range of A, where it is rn^ < 0, D~ and D^ protect the squared masses of 
charged sleptons from becoming negative. Nevertheless, the mass of the lightest state, ri is smaller 
than the experimental lower bound of 45 GeV coming from LEP I and it remains so also in the tiny 
region of A around 20 TeV, where tan/3 and /x are indicated by solid lines. 

When A increases, sneutrinos become heavier, whereas m~ tends initially to decrease. This is 
due to the fast increase of fj, tan/3 in the left -right entry of the r mass matrix, initially faster than 
the increase of the diagonal entries (see the A dependence of /X/M3 in Fig. |^. By identifying D^ 
and D~ (= D~), and neglecting third generation mass effects on the evolution of soft parameters, 
as well as Ar in the off-diagonal entry of the r mass matrix, one can cast the condition m^ < 

in the form: ml ml + D^{m,\ + m?^) — (/i tan/3 m,-)^ + D~ < 0. This quadratic equation in M2 
admits indeed two solutions which delimit an interval not too dissimilar from that indicated by the 
dotted portions of lines in Fig. ^ and determined without any of the above approximations. For 
further increases in A, the diagonal entries in the r mass matrix increase far more rapidly than the 
off-diagonal ones giving eventually a physically acceptable spectrum. 

We observe that choices of x closer to 1 (i.e. X ~ A), affect the super symmetric spectrum 
through the functions /(x), g{x) and the smaller size of the logarithm txQo- ^^ obtain lighter 
squarks and sleptons as well as lighter values of |?tt,|^ |: /x is, then, in general, ~ M2. All gaugino 
are heavier and therefore weak and strong contributions to B are also larger. In the range of A 
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Figure 4: Values of tan/3 -/Li obtained from tlic minimization of Vq- The dotted intervals indicate non- 
physical solutions. 



physically acceptable (A 1^ 50 GeV), the values of B are smaller than those obtained with our 
previous choice of x (see (|l8|) ): strong gauge couplings as well as top yukawa couplings are now 
larger at X, but not enough to compensate the decrease of txQo ^ue to smaller X's for fixed A's. 
(The values of Qo are also slightly smaller, but this change has negligible consequences.) Overall, 
the new tan/?'s do not differ appreciably from the values previously obtained, for A <; 50 GeV. 



3.2. Fully one— loop corrected Higgs potential 



We have previously argued about the suitability of our choice of decoupling scale Qoi which will 
be a posteriori justified by the results of the calculations presented in this subsection. It is obvious, 
however, that a change in this scale affects those parameters where the interplay between weak and 
strong sectors with different sensitivity to Qo has the largest effect. Particularly problematic is the 
parameter B. Since its low-energy value turns out to be rather small, changes of Qo may easily 
induce oscillations of B around zero, which in turn, require flips of sign in the parameter // with 
non-negligible consequences for the resulting phenomenology. Higher order corrections than those 
provided by the one-loop RGE as well as adjustments for the inadequacy of one unique decoupling 
scale Qo for the widely separated strong and weak sectors, become therefore mandatory. 

Given the smallness of the logarithm tq^x , finite corrections are in the MMM more important 
than in the MSSM, when compared to the leading logarithmic corrections supplied by the same 
order RGE. In principle, they should be also larger than corrections coming from higher order RGE. 
In practice, this will not be the case, for the parameter B. 

For a complete set of corrections to the parameters which determine the radiative breaking of 
SU{2) X U{1) we will consider the fully one-loop corrected effective potential Vi{Qo)^ defined as 
Vi{Qq) = Vo{Qo) + AT/i((5o) and we will include contributions to /S.Vi{Qq) coming from all sectors 
of the theory. Because of the wide gap existing between weak and strong mass parameters, this is 



crucial in order to warranty results stable under variation of the decoupling scale. 



Among other possible finite corrections, we include those to the b-quark mass, which, as known, 
can be rather large, when tan/3 has values as large as those obtained here. These corrections reduce 
the value of the coupling hb since hb{Mz) = {rub — 6mb)/vd = {rrib — tan/3 ^A)/(f cos /3) with A 
(> 0) given in [10,|ll|. Consequently, they affect the evolution ofrnj^^ (see ([A7D), decrease the low- 
energy value of B, therefore increasing tan/3. This, in turn, increases l/cos/3 and the size of 5mb- 
A retuning of all parameters finally produces the electroweak minimum for tan/3 's in general larger 
than those obtained when these corrections were not included. The values of ^ remain practically 
unchanged, since the parameters mj^ are unaffected by these corrections. For A = 100 TeV, in 
particular, it is /z ~ 397.5 GeV, as before; hb decreases from 0.798 to 0.786 and correspondingly 
At and Ab are slightly more negative. B adjusts to values indistinguishable from the previous one, 
but tan/3 increases to 52.1 from the initial value of 46.4. We neglect finite corrections to quark 
masses other than to nib, as well as threshold effects for supersymmetric parameters other than 
those induced by AVi{Qo)- 

We come now to analyze the effects of the inclusion of the one-loop corrections to the effective 
potential, which have the well known expression: 



AViiQo) 



1 



647r2 



E^amtiH) 



In 






(211 



where ma{H) is the field-dependent mass of the a -particle, n^, the corresponding number of 



degrees of freedom (negative for fermions): n~ = 6, 



-12, nr= 2, ni = -4, n~ 



-4, nH± = 2, 



nw 



6, n~Q 



-2, rijjo = 1, and nz = 3. The corrected potential Vi{Qq) yields minimization 



equations which retain the form in (|^), (^), provided the parameters Bfi, mjj , m'j^ are substituted 
by hatted parameters Bfi, fh'jj , mjj defined as mjj, = rnjj, + 5m\,, and B^ = B^ + 6{Bfj,): 



{ml^+fi^ + ^{g^+g^){vj-vl))2vd-(Bj,)2vu = 
{mj,,^+fi^-^{g^+g^){vj-vl))2vu-(Bji)2vd = 0. 



(22) 
(23) 



The shifts dmjj,, S{Bij,) are listed in Appendix B. We do not use different symbols to distinguish 
between Vu and Vd as obtained from the minimization of Vq and Vi. In what follows and in 
Appendix B we refer to the new minimum as to the Vi-minimum as opposite to the Vo-minimum 
obtained from (g). 

The corrections dm'jj, have been only partially included in @,P, but not those to B^, in 
principle very important and with strong impact on the value of tan/3 . In a generic GMSB, one 
can fix the value of tan/3 at will using the freedom in the high- and low-energy parameter B. It 
is precisely this lack of freedom which makes the MMM a model far more difficult to study. 

At the Vi-minimum, quark/squark and lepton/slepton contributions to 6{Bfi) are []: 
1 



327r2 



3 [At /u /i?) D(m~ ~ ) + 3 (a /i hi) D{m~^~) + [a, ^ h'^) D{m^^ 



those coming from charged gauge/gaugino, higgs/higgsino modes: 



1 



327r2 



2{M,^^g^)D{m~^~;, 



^ 2{Bii) + ^-v^siYi2(3\ D{mH±,G±] 



(24) 



(25) 



'Although squark mass matrices were evolved down in their 6x6 form, the approximation of 2 x 2 mass matrices 
for each generation of squarks and sleptons is used here 



where t;^ = f^ + v"^ and G^ is what would be the charged Goldstone boson at the Vo^nimimum. 
Those due to neutral gauge/gaugino, higgs/higgsino mode are: 



1 
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2 EF(m~o,Qo) ^1^ - ^^^ {{Bl.) + ^^^.^sin2/3) I)(m^o,ao) 



(26) 



with A^ the pseudoscalar Higgs, G^ the would-be-neutral Goldstone boson. The functions D{ma^b) 
and F{ma,Qo) are defined in Appendix B. Corrections to Bfj, coming from other sectors of the 
theory are identically zero. 



A comparison of (|24| ) with the RGE for B (see ([A4D ) shows how the quark/squark, lep- 
ton/slepton contributions to S{Bfj,) are corrections to the parameter B. They improve upon the 
arbitrariness of the scale Qq, linking it more realistically to the actual mass of the scalars virtually 
exchanged in the corresponding loop diagram. In the top/stop case, for example, in the limit of van- 
ishing left-right mixing terms, it is indeed, mj -^ m~ and m~ -^ ?ti-~, with [/, Q the two scalars in 

the first diagram of Fig. 3. Moreover, in the limit m~ — > m- = m~, it is D{mr ) — > 2log{'m^/QQ) 
and the term 3/(167r2)j4j/ij log(r?T^/Q§), with opposite sign to the corresponding one arising from 



(A4), would have the effect of trading Qq for m^. In its actual form, the top/stop contribution 



includes also finite corrections, of type c'-f at{X) a3{X)txQo i-e. with one logarithm less than the 
corresponding term provided by the RGE c^ at{X) as{X) (txQo)^- 

Similar considerations hold for chargino and neutralino corrections which take care of scale 
adjustments and finite pieces inclusions, of type ~ c'i"a{X), coming from the first diagram of 
Fig. |2|, to be compared to the terms cf a(X) txQo obtained from the one-loop RGE. 



The sign of the contributions (p4|)-(pq) depends crucially on the value of Qo with respect 
to the masses exchanged in the loop. Keeping the choice of Qq made at the tree-level, we have 
for A = 100 TeV, mj^ ~ 977 GeV, mj^ ~ 1099 GeV, and Qo ~ 993 GeV. These masses, roughly 
independent of further adjustments in tan/3 and fx, are shown in Fig. normalized to the gluino 
mass. Thus, Dlmr ) is positive; the same is true for Dim- ), whereas Dim- ) and Dim,- ) are 

negative. Quark/squark corrections are therefore small and negative {At and Ai, are both negative), 
those due to tau/stau are small and positive. Chargino corrections are large and negative. Although 
this cannot be explicitly seen by the formulas displayed in (|2g) and in Appendix B, negative are 
also the neutralino contributions to 6{Bfj,), although not as large as the chargino contributions. 
The remaining corrections coming from gauge and higgs boson modes are numerically smaller. 

Whereas the overall (5o~dependence of B^ and dm.'j^. is of order higher than the order of 
the calculation presented here, the relative size of the individual corrections changes for different 
choices of Qq. For smaller values of Qq, the quark/squark contribution tend to dominate and it 
may happen that approaching Mz, a negative // is needed to balance corrections too negative and 
therefore maintain Bfi positive {Bfi is indeed related to the fully one-loop corrected pseudoscalar 
Higgs mass). Since the bulk of supersymmetric particles is much heavier than Mz, we believe that 
definite conclusions cannot be drawn on the solutions obtained with such an "unnatural" choice of 
scale without including higher order corrections to the minimization conditions. 

As for the shifts drnjj., the two largest contributions are the top/stop correction to m|^ and 
the bottom/sbottom correction to m|^ . They are both negative and add to the already negative pa- 
rameters m^jj,. Contributions from first two generations squarks with different isospin, although not 
negligible in absolute value, almost completely cancel each other (see formulas in Appendix B). The 
next largest contributions are the chargino/charged-boson-sector and neutralino/neutral-boson 
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Figure 5: Solutions obtained for tan/?-/j, with a fully one-loop corrected neutral higgs potential (solid 
lines). The short-dashed lines correspond to the solutions obtained when only mj,-corrections are included; 
the dotted lines are the solutions in Fig. 0. See text for the long-dashed lines. 



sector, which contribute with different sign. Their algebraic sum, positive, is however still larger 
than the tau/stau contribution. The hatted parameters m|^ , m^ are now more negative than 
m\j and mhr . 

A retuning of /i and tan/3 is therefore needed to make the left-hand-side of (|2^ ) and (^) 
vanish. After a numerical iteration, solutions are found for values of tan/? and /i larger than 
those obtained at the Vo-minimum. Considerable is, in particular, the deviation for tan/3 : for 
A = 100 TeV, the new value is 59.0, the old, 46.4. This situation is not improved by a change of 
decouphng scale, which is symptomatic of some incompleteness in the set of corrections included. 

A close inspection of the two-loop evolution equation for B (see the relative equation in 
Appendix A) shows the presence of strong terms, not ^i-induced, which yield contributions of 
type do„ {as{X) at{X) txQo)^- The first of these terms, with n = 1, is of same size of the Ai- 
induced finite corrections. It is indeed originated by the same upper diagram in Fig. y, with the 
fermionic loop "open" (i.e. not shrank to reproduce At). We add this first term in our determination 
oi B; it has a positive sign and produces therefore an increase in the value of B and a decrease of 
tan/3. 



In complete analogy with the situation observed in the case of the one-loop RGE, the remaining 

c 
a^(X)txQo)" due to weak 



strong Aj-induced terms in the two-loop RGE yield corrections of type df „ {al{X) at{X)t'xQ^)'^ 



of same size than the corrections (ij^( 



and d^,,{as{X)at{X)tj,Q^r, 
terms. The first terms in these series (n = 1), are larger than the corresponding finite one-loop 
corrections, Cj^* as{X) at{X)txQo, c^" a{X), as expected. Their addition, therefore, would certainly 
require the inclusion of finite thresholds effects and corrections due to the two-loop RGE for all 
super symmetric parameters. No other anomalously large corrections are induced by terms in the 
two-loop RGE for m'jj and m'jj ||l 
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We show in Fig. ^ (solid lines) our estimate of /x, tan/3 as obtained following the prescription 
described above. We start our plots from A > 62 TeV []; below this value it is m^ < 0, and 
the calculation of one-loop corrections becomes impossible (see terms ln{m~/QQ) in the tau/stau 
corrections). We also neglect a very tiny region at small A which yields unacceptably low masses 
for Ti. For comparison, we give in this figure also the results obtained through minimization of Vq 
(dotted lines); those obtained when only finite corrections to rriij are added (short-dashed lines) 
and those coming from the minimization of Vi , with finite corrections to m^ included (long-dashed 
lines). In this last case tan/3 is large enough to keep mi negative up to A ~ 80 TeV. As already 
mentioned, corrections to rrib only do not produce deviations in the value of fi as obtained at 
the Vo-minimum (the dotted and short-dashed lines corresponding to ^ coincide). Similarly the 
inclusion of the non-resummed term in the two-loop RGE lowers tan/3 , but gives rise to changes 
in fj, indistinguishable on the left frame, and only barely visible in the right one, where the ratios 
n/M^ are given. 

The addition of all the above described corrections to the Higgs potential produce a much 
milder dependence of the tan/3-// solutions on the decoupling scale Qq. We give in the table below 
the values of tan/3 and // at the Vq- and the 14 -minimum for A = 100 TeV, and the two different 
decoupling scales Qo = 993 GeV and Qq = 331 GeV. 







Qo = 


= 993 GeV 




Qo = 


-- 331 GeV 


Vo-min 


tan/3 


= 46.4, 


// = 397.5 GeV 


tan/3 


= 36.1, 


^ = 481. 8 GeV 


Vi-min 


tan/3 


= 49.9, 


^ = 433.9 GeV 


tan/3 


= 47.3, 


fi = 446.0 GeV 



The large variation obtained when minimizing Vq can be traced back rather easily. The lowering 
of Qo induces an increase of the value of the strong gauge coupling in the interval {Qo,X}. This is 
due to the fact that the interval {Mz,Qo} where the faster SM evolution takes place (ftf*''^ = —7; 
63 = —3) is now smaller. The trilinear couplings At, A), are larger, and, as a consequence, also the 
values of B increase. Heavier squark spectra and larger top yukawa couplings induce also larger 
/Lf's; the increase of the term Bfj, produces smaller values of tan/3. The variation in the tan/3 -/i 
solution is much smaller when corrections to Vq are included. Values of Qq too close to Mz, are 
however still problematic. 

This shows that our choice of Qq is already close enough to the optimized scale Qs discussed 
in 1 13, 14], where the prediction for the two vacuum expectation values vi and V2 from Vq and Vi 
coincide. We can therefore rely on the RGE-improved tree-level scalar potential and our choice 
of Qo to verify whether dangerous minima breaking charge and/or colour are not also present and 
possibly deeper than the observed electroweak minimum. An analysis of this type when all one- 
loop corrections to the scalar potential are added would clearly be prohibitive. We do not find 



evidence for the existence of such ininima, at least using the analytical criteria given in |14.11]. 



4. Spectrum, Prediction, Constraints 

We display in Fig. ^ the spectrum obtained in the allowed range of A for our estimate of /i, 
tan/3 (solid lines of Fig. ^), and we give explicitly in the following table the values of gaugino, 
chargino. neutralino, and slepton masses (in GeV) at A = 100 TeV. 

"^This value of A is larger than the value of ~ 50 TeV obtained when minimizing Vb, since fi and tan/? are now 
larger 
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Ml 


M2 


M3 


m~+ 




Xi 


m~o 
X2 


m~o 
X3 


"^x2 


2^1,2 


m~ 

1^3 


m~ 




"^en 


"^e. 


144 


268 


757 


254 


456 


142 


255 


440 


455 


337 


333 


105 


361 


175 


347 



We follow here the conventions in [flSl (r?T,~+ > m~+) and we indicate by m~ (m~ ) the two 

IJ 1 ^ Xl X2 ^ •' BR '^ Sl' 

degenerate eigenvalues for right- and left-handed e and jl. Similarly degenerate are two of the 
three V states. For the same A, the up- and down-squark masses are given in Appendix C. Due 
to the rather large values of /x and tan/3 , the three heavy Higgs states (two neutral, one charged) 
range between 346 and 360 GeV; the lightest neutral state has mass roughly 125 GeV. 

The almost factor of two between the first two generations left-handed and right-handed 
masses, already present at high-scale (^) is clearly visible in the fourth frame of Fig. ^. The 
same splitting in the masses of the first two generations squarks is more modest since it is due 
to weak gaugino loops, whereas the bulk of the masses is produced by gluino loops ([T^). The 
corresponding ratios of the high-scale first two generations left- and right-handed squark masses 
(r77.^(X), m^AX)) over the low-energy gluino mass, not shown in this figure, are shifted downward, 

with respect to the solid lines in the two lower frames, by a factor ~ (1 — ZD/fts (for f{x)/g{x)'^ ~ 1) 



(see (13)). Their geometrical mean, Qq/M^, is explicitly indicated in both frames by the short- 
dashed lines. Tiny differences in the shape of the two solid lines in these two frames are due to 
isospin effects (12,^) as well as to the presence of different D-terms in the two cases. Those 



induced by intergenerational as well as chirality mixing terms in the up- and down-squark mass 
matrices are not visible in these figures. 

Right-handed as well as left-handed e and Ji can easily decay to the lightest neutralino Xi 
with mass ~ Mi . The decrease for increasing A of the ratios relative to the two heavy neutralinos 
(mainly neutral higgsino states) and the heavy chargino (mainly a charged higgsino state) in the 
first two frames of Fig. ^, can be simply explained by the milder growth of /x already observed in 
Fig. ^. Their masses actually increase and their states become increasingly more mixed with the 
gaugino states. 



Third generation yukawa effects in the evolution equations ( [A6| ) are responsible for the roughly 
10% decrease in the third generation squark masses. In addition, the two stop masses are only a little 

^1 



affected by the presence of left -right mixing terms in the up-squark mass matrix (~ Atrrit ^ Af| 



compared to the diagonal elements which are ^ M^). This effect is somewhat larger in the down- 
squark case. Differently than in the MSSM, the twelve squark states cluster quite closely around 
a common value of mass, ~ 1 TeV at A = 100 TeV. We show in Appendix C the diagonalization 
matrices for up- and down-squarks, at this same value of A, from where one can read off the 
composition of the relative squark mass eigenstates. 

Small is the effect induced by hr in the sneutrino sector, as the third frame of Fig. shows. On 
the contrary, the effects due to the presence of left-right mixing terms in the charged slepton mass 
matrix are rather large. The two states ri, 72, have masses on opposite sides of the left- and right- 
handed first two generations sleptons. The lighter of the two, ti, is indeed the lightest sparticle in 
the spectrum and exceeds the LEP I bound of 45 GeV only for A ^ 72 TeV. Once this constraint 
is imposed, the remaining spectrum is heavy enough to largely exceed any other experimental 
lower bound on masses. We obtain, in fact, rn~± ^ 180 GeV, m~o ^ 100 GeV, m~ ^ 240 GeV, 
mj^ k, 130 GeV, m~ ;^ 735 GeV, and mj ;^ 725 GeV. 

Thus, the main decay for Xi is the two body decay Xi ~^ t^i- It proceeds with full electroweak 
gauge strength reducing therefore the partial width for the decay Xi ~^ iG {G is here the gravitino) 
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Figure 6: Ratios of chargino, neutralino, and sfermion masses over weak and strong gaugino masses. For 
sleptons and squarks, tiie solid lines are relative to the degenerate first two generations left-handed (upper 
lines) and right-handed (lower lines) states. The dashed lines refer to the third generation eigenstates. The 
short-dashed lines in the two lower frames indicate the ratio of our choice for Qq (see text) over the gluino 
mass. This spectrum corresponds to the solutions (tan/3 , //) indicated in Fig. || by solid lines. 
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to practically negligible levels in this model. (For a theoretical discussion of this decay mode 
see [|l^.) The subsequent decay ri — > tG gives tt +missing energy as signature for the decaying 
Xi and a four r's signal for Xj* pair production in an e+e~ collider. Given the rather large values 
of m~o induced by the constraint m~ > 45 GeV, such a signal can be observed only at future 
e~^e~' colliders. The only signal LEP II is likely to detect is two r's +missing energy due to a pair 
production of ri. The two final r's are in this case more energetic than in the similar production 
mechanism in the MSSM, where ri decays into Xi^t ^^^ the neutralino has non-negligible mass. 

As already observed, an increase in x towards the value one, has the effect of decreasing scalar 
masses, while increasing gaugino masses. The parameter fi, now much closer to M2, gives chargino 
and neutralino states much more mixed: the curve relative to the heaviest neutralino in Fig. ^ is 
lowered to the values 1.3-1.1 in the shown range of A. The ratios of up- and down-squark over 
gluino masses are now between 0.9 and 1.0. The lightest particle of the spectrum is still ri, but 
also Br, Jill become lighter than the lightest neutralino. 

A severe constraint on this model may come from the 6 — > S7 test. (This was considered in |4,17| 
in the context of more general GMSB models than the MMM.) The estimate for the branching 
ratio Br{b — > 57) is obtained here using the one-loop supersymmetric boundary conditions at the 
electroweak scale given in Q and including the leading QCD corrections as in [^,18|. Experimental 



errors of some relevant low-energy variables as well as theoretical uncertainties in the QCD cor- 



rections are kept into account as in [18| and |19|. The main source of uncertainty comes from the 



ambiguity in the scale at which this process has to be evaluated: we let this scale vary between 

?TT,f,/2 to 2?TT,f,. 

Among the supersymmetric amplitudes, the collection of those relative to chargino exchange 
deserves a little attention in this model. We rewrite them in the form: 



A~. = C^5](-2)x^, [c\)k{Fi + -F2){x,k) + ^ Cl)k{F, + -F4)(x,,)| (27) 



2,6 

r2 /^2 



where the function Fi are given in p, the symbols Xij denote ratios of masses: x^^ = My^/mX , 

^fc 

Xjk = m~_^/m~ , and the constant C^ is C^ = Gi?e/(\/327r^). The coefficient Cj'!-^ collects the 
couplings: 

ribs _ ^jkb^jks* (ri'^bTTJks* TTJkb^jks*\ , TTJkbTTJks* /r,o\ 

of the "pure gaugino", "mixed gaugino-higgsino" and "pure higgsino" contributions. By using the 
definitions in 0] and the fact that sin/3 ~ 1, cos/? ~ l/tan/3 in this model, Cj**-^ can be re-expressed 



as: 



nbs r^ IT/ |2pfc3 ■p*fc2 Hhw^^ jr* f-pkS -pk3 ?>- ,-pk3 rk2 t^ \ , ^tw |^^ |2|Tfc3 |2 7<^ j^* (r,n.\ 
^Ijk - l^jll i~UL^UL-d—\^jlVj2[>'UL^UR^ts+i-UR^ULJ^tb)+ — \Vj2\ | i f/ijl ^tbJ<ts ^ i^^) 

where K is the Kobayashi-Maskawa matrix, U and V the matrices needed to diagonalize the 
chargino mass matrix. The 3x6 matrices Ff/^, and F^/^, juxtaposed, give the diagonalization matrix 
Dif of the up-squark mass matrix. Their numerical components as well as those of Tdl, ^dr, as 
obtained for A = 100 TeV, are given in Appendix C. Similarly, the coefficient C2% multiplied by 

m^Jrub is: 



Cl), - J^ tan/3 (-U,,V;, Vf.Tt'l + f-fU^.V;, V^An'^RKl) ■ (30) 
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Figure 7: Branching ratio for the decay 6 — > 57 obtained in the MMM (solid lines); the dashed lines indicate 
the SM's band and the dotted lines the band obtained in a Two-Higgs doublets model with the same mass 
for the charged higgs as obtained in the MMM. The horizontal solid lines delimit the experimentally allowed 
values. The results shown correspond to the solutions tan/3 -/i indicated in Fig. H by solid lines. 



The amplitude in ( p7D with the pure gaugino coupling (first term in (p9|)) is the supersymmetric 
counterpart of the SM amplitude: 



Asm = C^ KtbKl Sx*^ (-Fi + F2)[xt^) ; 



(31) 



those with pure higgsino couplings (the last terms in (|2g|) and (P0|)), the counterparts of the Higgs- 
mediated ampHtude: 



Ah- = C KlKtb xth jcot^ /? (^ ^1 + ^2) {xth) + (^ ^3 + i^4) {xth)^ 



(32) 



with obvious meaning of all ratios Xij. In particular, the amplitude with the last coupling in (pg|) 
corresponds to the Higgs contribution proportional to cot'^ j3, that with the last coupling in (RO) to 
the tan/3 -independent one. Both couplings have column-index equal to 2 in the diagonalization ma- 
trices U and V , which, indeed, selects the higgsino component of the chargino exchanged in the loop 
(see [|lq ], whose notation we follow here). We remark that in this model the chargino mass matrix, 
X ~ ((M2, -v/2Mvi/), (0,/x)), has all positive entries. It is, therefore, detX > 0. The two diago- 
nalization matrices C/, V have then both the structure [{cos(j)uy ,sm. (I)uy),{— s\n (puy, cos (puy)) 



with angles 
,2 



(j)v determined by tan( 



[m~ 



XI 



(m~ — M2) / {\f2Mw M2) , both > 0. (We remind that in our convention 
chargino.) 



{M'i + 2M^))/{^/2Mwf^) and tan (py ~ 

is the heavier 



m 



XI 



The two amplitudes with couplings ( |30| ) are, in general, those responsible for the growth of 
the chargino contribution for increasing tan/3 and should, in principle, be the most relevant ones 
when tan/3 is large. The situation is, however, slightly different in this model. As it can be seen 
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from Appendix C, the tan/3 -dependent gaugino-higgsino amplitude, 

gets contributions from the mainly left-handed stop {k = 2) and the mainly left-handed scharm 
[k = 6), which have non-negligible mixings with cl and t^, respectively. (The index k is ordered 
according to increasing values of the squark masses.) Sizes and signs of these mixings are determined 
by the RGE which generate off-diagonal terms in the left-left sector of the up-squark mass matrix, 
and the requirement of orthogonality for D. For our representative value A = 100 TeV, ir^^^rj/^l ~ 
I-'^cFl-'^c/lI ~ 0.02. The two products have opposite sign and differ only at the 10~^ level; the 
presence of the functions (F3 + 2/3F4)(xj_fc) ~ C(l), however, gives rise to a flavour violation 
still at the 10~^ level. An additional cancellation is due to the gaugino mixing. After summation 
in /c, the sign of the two terms corresponding to the two different j's is different, as it appears 
evident from the previous discussion on the matrices U and V . (This result, obviously does not 
depend on the freedom used in fixing relative phases in the elements of U and V (see fl^).) Again, 
for A = 100 TeV, it is: Ui2V{. ~ -0.36, U22V^. ~ +0.20, (m~ = 456 GeV, m~ = 254 GeV). 
Furthermore, the ratio x^fc is a strong factor of suppression, x^k ~ 10~^, reducing the "gaugino- 
higgsino" mixing amplitude ( |33D to be more than two orders of magnitude smaller than Asm and 
Ajj- and roughly one order of magnitude larger than the pure "gaugino" contribution. Similar 
cancellations appear also in the other gaugino-higgsino mixing amplitude (with coupling in the 



second term of (29)), which gives rise to the smallest amplitude in the collection (| 



The tan/3 -dependent pure "higgsino" amplitude gets non-negligible contributions from the 
two stop states (A; = 1 and k = 2). The product of the two elements of Tjjl, Tur, with their 
opposite sign for the two /c's, leads to a left -right mixing < 10~^, overkilling the enhancement due 
to tan/3 . Numerically, once all summations on k and j are performed, this amplitude turns out to 
be smaller than the tan/? -independent pure "higgsino" amplitude: 

"^HH = C^Y.(-2)x^k^\V,2?\T''^R?Kt,Kl(F, + lF2){x,k). (34) 

3,k 

where it is (F^^)^ ~ 1, V'22 — ^ ^'^'^ ^tw/2 — 1 and where the flavour violation has the same weight 
as in the SM. The amplitude (|3^ is the largest one in the collection (|27|). Its sign is opposite to 
that of Asm and ^//-, but due to the suppression factor x^k-, it produces only a tiny cancellation 
oiAH-. 

Also the gluino amplitude: 

^5 = ^" E IT ^""^ (— ) [^DL ^DLX*^' F^ixgk) - ^ T% T*,^l F,{x,k)\ (35) 



is known to increase for increasing values of tan/3 [20|. Of the two terms in (p5|), the larger 
contribution comes from the second one, with exchange of the mainly right-handed sbottom (A; = 1), 
the mainly left-handed sbottom (k = 4) and the mainly left-handed s-strange (A: = 6). Left -right 
terms in the down-squark mass matrix play a larger role than in the up-squark mass matrix 
and the corresponding mass eigenstates are more mixed states. The simultaneous left-right and 
flavour transitions are, however, still ^ 10~^. The enhancement factors mg/mi,(mb) and Us/a^, 
together with the fact that down- and up-squarks have very similar masses, makes this amplitude 
comparable to the tan/3 -dependent pure "higgsino" one. In our representative case A = 100 TeV, 
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the SM, Higgs, chargino and gluino amplitudes are respectively: —1.1 x 10^^, —5.7 x 10^^'^, 7.9 x 
10~^^, and 1.3 x 10^^^. The charged Higgs mass, for this value of A, is 360 GeV. The neutralino 
amplitude is completely negligible. 

The values of Br{b -^ sj) predicted by this model are in the region delimited by solid lines in 
Fig. 1^ The MMM's branching ratio deviates little from that obtained in a Two-Higgs-Doublet 
Model (band delimited by dotted lines) with mjj± as predicted by the MMM. For reference we 
report also the SM prediction (region within the dashed lines) and the experimentally allowed 
value (region within the two horizontal solid lines) pl|]. 

Since the requirement of positivity for the the ri-mass has already selected relatively large 
values of A, A ^ 62 TeV, and the LEP I constraint m~ ^ 45 GeV imposes A ^ 72 TeV, no further 
exclusion comes, at the moment, from the measurement of this decay. Improved experimental 
results will require improvements in the calculation of the supersymmetric branching ratio before 
any significant conclusion can be drawn. A significant shift of x towards one, however, would push 
already now the MMM's band outside the experimentally allowed range up to A ~ 100 TeV. For 
this value of A, it is: mfj± = 300 GeV and the SM, Higgs, chargino and gluino amplitudes are: 
-1.1 X 10"^, -6.7 X 10~i°, 4.9 X 10"^2 and 1.5 x 10"^^ 



5. Conclusions 

In this paper, we studied the MMM in the approximation of a messenger scale of the same order 
than the supersymmetry-violating messenger scalar mass, and gave our estimate for the solutions 
tan/?-/x enforced by the mechanism of radiative breaking of SU{2)l x [/(l)y. 

Differently than generic GMSB models, which benefit from the presence of more free parame- 
ters at the messenger scale, the MMM requires a careful handling of corrections to Bfj,, mjj and 
mjj , which determine the electroweak breaking. To this aim, we considered the one-loop corrected 
effective potential with contributions coming from all massive modes in the model: quark, squark, 
lepton, slepton, chargino, neutralino, gauge and higgs modes. Finally we included additional loga- 
rithmic two-loop corrections to the parameter B of same type and size than those induced by the 
one-loop effective potential. All of these corrections turn out to be important to obtain solutions 
tan/3-// stable under variation of the decoupling scale Qo around a typical squark mass. 

Among the predictions obtained by studying the MMM's mass spectrum, the most interesting 
is that Ti is the lightest sparticle. Indeed, the requirement of positivity of m^ as well as that 
coming from the LEPI lower bounds on supersymmetric masses, exclude already values of A up 
to ~ 70 TeV. This prediction has other important consequences. The MMM, in fact, a) cannot 
accommodate the e"'"e~77 CDF event |^]; b) may be detected in e"'"e~ collisions through two r's 
+missing energy or four r's +missing energy signals. 

As expected, flavour violation effects as well as effects due to left-right mixing in the squark 
sector are small in this model. The sparticle contributions to 6 — > 57 , opposite to the W and H 
contributions are, therefore, not very significant. The overall MMM's prediction closely resembles 
that of a Two Higgs Doublet Model with same H^ mass. When m~ is larger than 45 GeV, we 
obtain rates consistent with the present experimental measurement. 
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Note added 

After completion of this work we became aware of the existence of the paper hep-ph/9701341 by 
D. Dicus, B. Dutta, and S. Nandi which studies one of the experimental signals listed in this paper, 
and of the content of hep-ph/9612464 by R. Rattazzi and U. Sarid with the same subject of the 
present paper. The solutions tan/3 -/x seem in qualitative agreement with ours. Differences appear, 
however, in the conclusions reached for the implication of the measurement of 6 ^ 57 . 
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Appendix A. 

We list the one— loop RGE to which we refer in the text. For simpUcity we give the approximate 
form vahd in the hmit of Kobayashi-Maskawa matrix i^T ~ H. The one-loop equations actually 
used in our analysis, with a realistic K, are given in y]. 

— gauge couplings and gaugino masses: 

Ui = -h, af Mi = -bi aiMi (Al) 

The dot indicates derivative with respect to txQ = 2 In {X/Q) and the tilde a division over 47r. The 
super symmetric /^-function coefficients bi are (—3, 1,33/5). The corresponding SM ones, bf^^, are 
(-7,-19/6,41/10). 

— yuka^va couplings and soft trilinear couplings: 

at = iGY)^ - 65t - 5fe A = - {GA)^ - 6Atat - Abon, 

ah = {GY)^ - Gofo -at-Ur A^ = - {GA)^ - GAbOb - Atat - ArOr 

5^ = (GY)^ - Aur - Sab A^ = - (GA)^ - AA-rUr - 3^5,, (A2) 

where at = hf/Air and after defining G^ = (13/36), (7/36), (3/4) for i = t,b, r, GY and GA are: 

(Cn = 4{C353 + C252 + |cf ai} 

5 

iGA)i = 4{C353M3 + C252M2 + |cf aiMi} (A3) 

5 

— bilinear coupling: 

B = -A\C2a2M2 + -(-) 5iMi} - (^^5^ + 3^5;, + 2,Atat) (A4) 

— first tvi^o generations squark and slepton masses 

K)ii,22 = 4 |c3a3M| + C2a2M| + ^Y^a^MlX (A5) 

— third generation sfermion masses: 

(my ={rnV) - at {S S\ - ab {S S\ 
Q 33 Q 11,22 

(ml) = (m^) -2at{SS), im'^-) = im^ -2ab(SS). 

(ml) = (ml) - ar iSS)^ (ml) = (ml) - 2 5^ (SS)^ . (A6) 

^ L^33 ^ L^l,22 r\Jr \ ^J ^^ \ £;/;^^ 22 ^ \ Jt \ J 

— soft Higgs-potential parameters: 

WJ = (4)n,22 " ^"^ (^^^^ " "- ^^^)- ' ^""'^"^ = (4)ll,22 " ^"* ^^^)* (^^) 

where 

{SS)^ = {ml + ml)33 + ml^+A^t 

(55), = iml + ml)33 + mjj^+Al (A8) 
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H parameter: 



(/^^) 



Q 1 9 

U C2OL2 + -i-j aij - (3at + 3q6 + Qt 



t(C72 



M 



(A9) 



The low-energy parameter ^ is obtained from the minimization condition of the scalar Higgs 
potential. This equation may help to trace back the high-energy value of this parameter. 

We report also the additional two— loop terms in the RGE for B to which we refer in the text 



(see HI): 



B 



80 1 2 207 

— (C2 a2f M2 + — {C2 02) (5i) (Ml + M2) + -- {aif Mi) 
00 25 



+12/ 


C3 a^M^ + - 

it. 


-aiMi 


at + 


C3 asMs - 


^aiMi 


Oib + — «iMia^ 


A 


C3 as + — ai 


Atat + 


CaSs-^ai 


Aat + — "l^rOr^ 


+ 6 


Aral + {At + 


Ab) atab + 3Ataf + SA^al 







(AID) 



Appendix B. 

It is easy to see that, if we neglect phases other than that in the CKM matrix, the scalar 
potential V (here indifferently Vq or the fully one-loop corrected potential Vi) depend on the real 
components (j)i,ipi {i = u, d) of the two neutral higgses H^ and H^ {H^ = (pd+iipd , H^ = (l)u+iipu) | 
through the bilinear operators 



erf 



I u' 



0|2 



\H. 



0|2 



£3 



HM + h.c. 



(Bl) 



(If extra phases are allowed, one would have to consider slightly more general operators. The 
following results would remain unchanged up to complex conjugates of soft parameters and/or /i). 

The derivatives of V with respect to </>j, which determine the minimum conditions, can then 
be expressed in terms of derivatives with respect to e^, £« and £3 as follows: 



{d^^V) = 2<A„(a,„y) + 20^9,3 y). 



(B2) 



The derivatives de^V are closely related to the entries of the pseudoscalar Higgs mass matrix, if one 
imposes that (V'd) = {i^u) = 0: d'^V/di^l = 2d,^V; d'^V/di^dd^Ju = -2de^V; d'^V/d^l = 2d,J^ . 
Moreover, at the electroweak minimum, i.e. when d^^V = dffy^V = for {(l)d) = Vd, {(pu) = Vu 
(tan/3 = Vu/vd), it is: 



d.,V 



-d,,Vtan(3; d,^V = -d,,V cot P . 



(B3) 



In other words, the correct minimization condition guarantees the existence of one Goldstone boson 
mode, corresponding to the eigenstate (cos/3, — sin/3). The massive mode (sin/3,cos/3) has mass 
(2/ sin 2/5) (— (9^3 y). In this sense —d^^V identifies {B^) when 1/ is Vq and what will be called Bfi 
in the case of Vi . 



^The physical neutral Higgs field iif°, /i", A", properly normalized, are expressed in terms of H^ and H'^ as in 
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The derivatives 5^^^, 3^^, and 9^3 of the RGE-improved tree-level potential, are: 



ml^+^l^-{g^+gl)e^/4: 



(B4) 



with e± = €(1 i Gu- Those of the one-loop corrected potential Vi maintain the same form when 
expressed in terms of "hatted" variables injj , m^ , and Bfi, where m|^. = mjj, + drnj^., and 

Bfi = B^+5{B^). The shifts are due to the one-loop corrections coming from the following sectors: 
quark/squark, lepton/slepton, chargino/charged gauge and Higgs mode, neutralino/neutral gauge 
and Higgs mode. They can be evaluated generalizing methods used in p4| , p5[ . We list them in the 
following. 

top stop contribution 



5m 



HdHA 



327r2 
3 



S[mr_ r. ) tt-^ 1- lJ{m- 



tl,t2' 



tiM' 



^h,^,f + il+EYlc^^2 



5m\, I , r — ?r 



{g'+gi) 



.2 I „2 



^2 {g'+gy) 



^K,J.) [h't-'^^^^] +Dim~:rJ [ihtAtr-^^^-^C^Af\-2F{mt,Qo) 



^(Bf^Al 



327r2 



D 



(""tlJa)!^*/^^?) 



(B5) 



where C~ is C~ = (1- fsin^^H^) and Al is defined as A^ = (m~ - m~)33 + e_ C~(g2+5'y)/4. 
The functions S{mj j ), D{mj j ), where the dependence on the scale Qo has been suppressed for 
compactness, can be expressed in terms of F{m, Qq) = 2m? \\og{m'^ / Q'q) — 1] as: 






(B6) 



The approximation of 2 x 2 mass matrices for each generation of sfermions is used here, in spite of the 
fact that the building blocks for these matrices are determined without neglecting intergenerational 
mixing terms. All the corrections listed in this appendix are clearly not sensitive to these tiny effects. 



bottom sbottom contribution 



5m, 

5m, 

S{Bii)\ 



2 


b,b 


3 


Hd 


32^2 


2 


fe,I ~ 


3 


Hu 


327r2 


',A 




3 



-.2 I „2 



'5("^, 



bi,b2' 



"^K,!, 



\lJ_l±^\+Dim~jj(ih,A,f- 



ig^+gy) , n^_ ^hu .a2 , (5'+5y) ^ ^2 



±^OrAA-2Fim,,Qo) 



b'-^b 



+ D{m~^~,U{hbf^r + 



CrA4 



b^b 327r2 



^K„62)(^^^^b) 



(B7) 



with C~ = (1-| sin2 6w) and A^ now defined as Z^ = (m?- - m~)33 - e_ C~(fi(2+^2^)/4. 

sneutrino contribution 

For each sneutrino species, we have: 



5m 



1 



H. 



327r2 



F{m~,Qo) 



ia'+gl) 
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5(5^)1- 



1 



327r2 
0. 



-F(m~, Qo 



(5^ +5?-) 



(B8) 



tau stau contribution 

It can be obtained from the bottom-sbottom contribution after suppression of the colour factor 3, 
with obvious replacements, i.e. h^ — > hr, nib -^ rrir, mr -^ ttt.- , C~ -^ C~ = (1 — 4sin2 Oy/), and 
A| - Af ^ (mi - m|)33 - e- C~(g2+g2.)/4. 

first t^vo generations quark/squark, lepton/slepton contribution 

Since the approximation of massless first-two-generations quarks and leptons is used in this anal- 



ysis, this contribution can be obtained from (B5,B7) and the tau-stau contribution putting to zero 
yukawa couplings, quark and lepton masses, and substituting third generation sfermion masses with 
the corresponding first and second generation ones. This contribution is rather small: there are no 
corrections coming from this sector to S/x, whereas the corrections to m?^ and m?^ from the two 
component of a SU{2)l doublet supermultiplet, although singularly not small, tend to cancel each 
other without the diversifying effect of the yukawa couplings. 

chargino/charged boson contribution 

The spin 1/2 modes, charginos, contribute as: 

1 



5mH, 



d'X^ 



^^hJ^± 



6{B^^)\- 



327r2^ 

1 2 

327r2^ 



g'e- 



1 



■9 



whereas the charged gauge boson mode gives: 



-^(-x.x.)-^K.x.)(^l + ^'-5'^ 
Z)(m~ ~)(2M2/x) 



(B9) 



(5m 



hJw^ 



5m 



hJw± 



327r2 

3 

32^ 



5{Bii)y± = 0, 



F{Mw,Qo 



F{Mw,Qo)y 



(BIO) 



with Mw defined as M\y = e^g^ /2. There are also two charged Higgs modes H^, G^ with "masses" 
given by 



m 



H±,G± 



1)2,2,5^ , 



^L- 



2 , y 



+ 2(i?/x) + |-e3 



with fMjj, = m'jj, + /u2. At the Vq- minimum, where the relations f^jj + {g'^+gY){v'^—v'^)/4: = 
(i?//)tan/3 and /x^-^ — ((7^+5y)(f^ — fu)/4 = (i?//)cot/5 hold, G^ is the charged Goldstone boson, 
H^ the usual charged Higgs, with mass m^± = //|^ +A*fi' +M^^. Both modes, however, contribute 
to the determination of the Vi -minimum (where the charged Goldstone boson, as well as the neutral 
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one mentioned below, are not massless anymore). They shift the variables Tiijj , m'j^ and Bfi by 
the quantities: 



^^hJh±,g± 



1 



6m 



hJh±,g± 



KBf^)\H±,G± 



327r2 
1 

1 

32^ 



4 



j'i 



S{mH±,G±) X + D{mH±,G±) V ( "^L ~ ^Hu + ^e- 



4 

^'2 



2 

-,'2 



5'("^H±,G±) -T - D{mH±G±) -r- ^Ha " '^H^ + ^e_ 



-^("^//±,G±)^(2(^/^) + y^3 



(Bll) 



neutralino/neutral boson contribution 

In this sector, there are: i) one gauge boson mode, with "mass" Mz = e+(5(2 + (/y)/2, which 
contribute as: 

ig'Wy)' 



^^hjz 



^mljz 



327r2 
3 



F{Mz,Qo) 
F{Mz,Qo) 



4 

ia'+g'y) 



327r2 
d{Bfi)\z = 0; 
ii) two neutral pseudoscalar Higgs modes A^,G^ with "masses" 



(B12) 



^AO,G0 = 2 ^/^L+^k± 



/^L-^k+^^^-l +mBf^)y 



At the Vo-minimum G^ is massless, A^ is the usual pseudoscalar Higgs with mass 
mjj + 2/^2 . They contribute as: 



^Ao = ^H, + 



^hJao,go 



1 (g'+gi^) 



327r2 8 



5m fT 

,2 , 1 (g'+g'y) 



DimAO^Go) mjj^ - mjj^ + 



(g'+gi^) 



5 m 



HjAO,GO 



327r2 8 
SiBl^)\AO,GO = 0; 
iii) two neutral scalar Higgs modes with "masses" 



Dim^ono) mjj - mjj + 



{g'+gi) 



(B13) 



2 1 J 2 ^ 2 , (g'+gy) , 



^^H,-^^HA{g'+gyy^] + 2(i?Ai)+ 



(5^5^) 



es 



which, at the Vo-minimum, reduce to the conventional neutral Higgs masses m?^(, ^o = {m|^ ~^^^'H ~^ 
Ml ± [(/u?^^+/u^„ -M|)2 cos2 2/3 + (^1^^+^^^ + M|)2 sin2 2/3]i/2}/2. They produce the shifts: 



<^"i//dko,/iO 



HBf^)\HOhO 



1 


ig' 


+gi) 


327r2 
1 


ig' 


4 

+gl) 


327r2 
1 


{g' 


4 
+5?-) 






327r2 



-D{mHo^h<->) (BfJ-) + 



(5^+5^) 



£3 



(B14) 
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iv) four spin 1/2 modes, the neutralinos. We find their contribution following [26,2^,| 

1 



•^"^Llp 



'^"^klxO 



S{Bf^)\~o 



32^2 

1 

32^ 

1 

32^ 



-2^F(mp,Qo)5,,m|o 

1 = 1 
4 

-2 5]F(m~o,Qo)5,„m|o 

i=l 
4 

+2 V F(m~o, Qo) 5^3 m~o 



4 = 1 



where in turn, the derivatives of the masses are 



de.m^o 



Ai A-i A,i 

cq + 2c4 mlg + 3c2 mig + 4co m^-o 

Xi Xj Xj 



ie^ 



^d^'^u^^S) 



(B15) 



(B16) 



The coefficients Cj are the coefficients of the characteristic equation com o + C2 ?7i o + C4m o + 

Xi Xi Xi 

cqw?q + cg = 0, for the matrix M.~oMX with all Higgs dependences reinstated back in M.~o- 
Xi X X ^ 

(The fact that at the true minimum charged Higgses do not develop vacuum expectation values is 
already imposed, allowing therefore the chargino and neutralino mass matrices to decouple.) They 
are : 



Co 

C2 
C4 

C6 

C8 



+1 

-2^?-Ml-Ml-{g^+gl.)e^ 

+/ + 2^2(^2 + Ml) + MlMJ + [^?{g^+gl) + {MJg'^ + Mlg^)]e+ 

-/x(Mi<7'2 + M2<7')e3 + {g^+glfe+V^ 



-^i^{Ml + M^) - 2^i^MlMl + /i[/i2(Mi5'2 + Ms*?^) + MiM2{M2g'^ + Mig'^)]e3 

-f^Hg^+gyfedeu - fi^M^g'^ + Mfg^)e+ - {M2g'^ + Mig^fe+^/A 

+f,^[^i^MfMi - ^MiM2(M25'2 + Mig^)e3 + {M2g'^ + Mig^feaeu] • (B17) 



Appendix C. 

We collect here the numerical entries of the up- and down-squark mass matrices, their eigen- 
values and eigenvectors, as obtained at A = 100 TeV. For better readability, we write here "linear" 
mass matrices: their elements are the square root of the entries in the real mass matrix; the re- 
ported signs are those of the corresponding squared elements. The up-squark mass matrix, split 
in the four left -left, left-right, right -left and right -right submatrices is: 



Mu 



/ 1099.432 


0.000 


-25.320 


0.000 


0.000 


-3.879 \ 


0.000 


1099.432 


53.712 


0.000 


0.000 


8.228 


-25.320 


53.712 


1044.930 


0.000 


0.000 


38.766 


0.000 


0.000 


0.000 


1054.342 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1054.342 


0.000 


^ -3.879 


8.228 


-38.766 


0.000 


0.000 


976.915 / 



(CI) 



As in 1^], we work in a quark basis in which the transition from current to mass eigenstates is 
obtained through the rotation of the left-handed up quarks only. The rotation matrix is the 
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Kobayashi-Maskawa matrix. The effect of this matrix on both up- and down-squark mass matrices 
is kept up to 0{X^), where A is the Cabibbo angle. 

The diagonahzation matrix for Mu is: 

/ .0001 

.0053 

.0000 

.0000 

-.9762 

V -.2169 



Du 



0004 


.0108 


0238 


.9996 


0000 


.0000 


0000 


.0000 


2169 


.0000 


9759 


.0244 



0000 


.0000 


.9999 \ 


k = 1 


f 976.907 \ 


0000 


.0000 


-.0108 


2 


1044.902 


0369 


.9993 


.0000 


3 


1054.342 


9993 


.0369 


.0000 


4 


1054.342 


0000 


.0000 


.0000 


5 


1099.432 


0000 


.0000 


.0001 / 


6 


^ 1099.466 / 



(C2) 



Each row k represents the decomposition of the ^-eigenvector over the initial basis {ul, CL,tL, ur, cr, Ir) . 
The 6x3 matrices on the left and right of the vertical bar are, then, nothing else but the ma- 
trices TuL and Tur used in the text. The last column of numbers in ( |C2| ) gives for each k, the 
corresponding eigenvalue in GeV. 

The mass matrix for the down-squark sector looks like: 



/ 1102.331 


0.000 


-30.548 


0.000 


0.000 


0.000 \ 


0.000 


1102.331 


64.803 


0.000 


0.000 


0.000 


-30.548 


64.803 


1033.967 


0.000 


0.000 


-256.273 


0.000 


0.000 


0.000 


1050.701 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1050.701 


0.000 


\ 0.000 


0.000 


-256.273 


0.000 


0.000 


999.164 / 



Md 



and the corresponding diagonahzation matrix is: 

/ .0018 -.0081 .5101 

.0000 .0000 .0000 

.0000 .0000 .0000 

-.0073 .0329 -.8595 

-.9762 -.2169 .0000 

V .2168 -.9756 -.0332 



Dd 



.0000 


.0000 


.8601 \ 


k = 1 


/ 979.331 \ 


.9998 


.0211 


.0000 


2 


1050.701 


.0211 


.9998 


.0000 


3 


1050.701 


.0000 


.0000 


.5100 


4 


1052.701 


.0000 


.0000 


.0000 


5 


1102.331 


.0000 


.0000 


.0101 / 


6 


^ 1102.397 / 



(C3) 



(C4) 



where the 6x3 matrices on the left and right of the vertical bar are Tdl and Tj:)r. The differences 
in the elements (1,1), (2,2), (4,4), (5,5) in Mu and Mjj are due to the different D-terms present in 
the two sectors. The off-diagonal elements (1,3), (2,3) in M^ are induced by the non-diagonal up 
quark yukawa couplings during the evolution from X to Qq. The same elements in Mu sum also 
the supersymmetry non-violating terms induced by up quark yukawa couplings. 
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